Elastic scattering loss of atoms from colliding Bose-Einstein condensate wavepackets 
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Bragg diffraction of atoms by light waves has been used to create high momentum components 
in a Bose-Einstein condensate. Collisions between atoms from two distinct momentum wavepackets 
cause elastic scattering that can remove a significant fraction of atoms from the wavepackets and 
cause the formation of a spherical shell of scattered atoms. We develop a slowly varying envelope 
technique that includes the effects of this loss on the condensate dynamics described by the Gross- 
Pitaevski equation. Three-dimensional numerical calculations are presented for two experimental 
situations: passage of a moving daughter condensate through a non-moving parent condensate, and 
four-wave mixing of matter waves. 



A light-induced potential applied to a Bose-Einstein 
condensate (BEC) can be used to make high momentum 
daughter BEC wavepackets which propagate through the 
parent condensate 0-0]. High momentum means very 
large in relation to the mean momentum in the par- 
ent wavepacket and the momentum mv s where v s is the 
sound velocity in the parent. Such techniques have been 
used to make an atom laser JW, to study the coherence 
properties of condensates |3],|5|,|6| , and to study nonlinear 
four- wave mixing (4WM) of coherent matter waves j^D . 
As explained in this Letter, elastic scattering between 
condensate atoms from different momentum wavepackets 
can remove copious numbers of atoms from these mov- 
ing wavepackets. Recently, profuse elastic scattering of 
atoms between daughter and parent BEC wavepackets 
has been observed at MIT || . Such losses will be an im- 
portant consideration for atom optics applications. Fig. |l| 
schematically shows wavepackets in momentum space 
where the high momentum wavepacket with central mo- 
mentum P2 = 2Kk. p h was produced by optically-induced 
Bragg scattering from the Pi = initial wavepacket. 
Here hk p h = h/X p h is the photon momentum for light 
with wavelength X p h- The spherical shell in Fig. [I] (ex- 
cluding the condensate wavepackets) results from elastic 
scattering between atoms from the Pi wavepacket and 
atoms from the P2 wavepacket. The elastically scattered 
atoms in the spherical shell can neither be described as 
part of the mean-field of the BEC, nor can the formation 
pot or evolution of the spherical shell be modeled using 
the usual Gross-Pitaevskii equation (GPE) pit , Eq. (Q). 
In what follows we use the term "elastic scattering" to 
mean only those non-forward elastic scattering processes 
not accounted for within the GPE. 

Here we provide a simple means of describing the loss of 
atoms from the condensate wavepackets due to the elastic 
scattering mechanism. This is made possible by (a) using 
the appropriate momentum dependence of the nonlinear 
coupling constant in the GPE [ ^2|JT^ ] and (b) using a 
newly developed slowly-varying-envelope approximation 
(SVEA) for the condensate wavefunction in systems with 



both slow and fast momentum components JL4J . Since 
the SVEA treats each distinct momentum wavepacket 
separately, we can incorporate the correct momentum de- 
pendence in the nonlinear coupling constants. Thus, we 
can treat the effect of elastic scattering losses on the con- 
densate dynamics using the SVEA version of the GPE, 
even in single spin component systems. We first outline 
the theory for describing elastic scattering loss and then 
present two examples, one applied to output coupling of 
atom laser pulses from a BEC source, and the other to 
4WM. A more complete discussion of the theory and fur- 
ther applications will be presented elsewhere p5[. 
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FIG. 1. (a) Wavepackets in momentum space with lab- 
oratory frame central momenta Pi = and P2 = 27ikph 
undergoing elastic scattering to produce a spherical shell of 
elastically scattered particles. In the CM frame moving with 
velocity v cm = 7tk p ;,/m, momenta are shifted by — h\i p h- 

The GPE for a single spin component BEC at zero 
temperature can be written as [ p~T]| 

ih^ = (T x + V(x,t) + U N T \V\ 2 )y , (1) 
at 

where T x = V 2 is the kinetic energy operator, V(n, t) 
is the external trapping potential imposed on the atoms, 
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Nt is the total number of atoms in the condensate, and W 
is normalized to unity. Although the coupling constant is 
usually expressed as Uq = 47r ^ t a ° , where ao is the s-wave 
scattering length and m is the atomic mass, it is more 
correct at zero temperature to express Uo in terms of the 
many-body T-matrix, which is often well-approximated 

by the 2-body T-matrix g||§: U (k) = - 4 "^ (fc) . 
Here T — S is the unitary iS-matrix, and Tik is the 
relative momentum of the colliding atoms. If we assume 
no inelastic scattering and expand T(k)/k in powers of k 



using S(k) 
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In a normal condensate, the lead term in Eq. (Q) gives the 
usual GPE. However, if two wavepackets with very differ- 
ent momenta interact (cross-energy terms in Eqs. (|^-0)), 
kao need not be negligible (kao = 0.06 for the Na con- 
densate example given below), and the second term in 
Eq. (||) must be taken into account. This term generates 
the elastic scattering loss term in the GPE. If we apply 
the optical theorem to the forward scattering amplitude, 
it is this second term that is responsible for scattering 
out of the forward direction. Suitable generalizations of 
Eq. (||) allow inclusion of inelastic collision losses and 
treatment of multiple spin-component condensates. 

The application of an optical standing wave pulse 
diffracts a fraction of the initial condensate into a high 
momentum component To a very good approxi- 

mation, the wavefunction immediately after application 
of a set of optical pulses is given by the superposition of 
wavepackets, 



J 

4'(x, t = 0) = ^>(x) 2^ a i ex P(*kj • x ) ; 



(3) 



where J is the number of distinct momentum wavepack- 
ets present. Here is the initial wavefunction of 
the parent condensate before application of the optical 
pulses; it is the solution to the GPE with a harmonic po- 
tential centered around x = 0. We assume the momen- 
tum differences ?i | — kj- 1 to be much larger than both the 
momentum spread in the initial parent BEC wavepacket 
and the momentum mv s associated with the speed of 
sound v s in the BEC (hence, superfluid suppression of 
collisions |j] does not occur). Since different wavepack- 
ets do not overlap in momentum space, J^f =1 \o,j\' 2 = 1- 
The SVEA is made by writing the wavefunction as 



much slower scales than the phases. Consequently, full 
three-dimensional (3D) calculations of the envelopes is 
numerically tractable, as we describe in more detail else- 
where 1^,0. In the first example we consider below, 
we take only two components, i.e., j — 1,2 and the ini- 
tial condition is $j(x,i = 0) = a,-i/>(x). In the 4WM 
process, j = 1,...,4, with $j(x, t = 0) = a,jip(x) for 
j = 1,2,3, and the j = 4 envelope is initially unpop- 
ulated, $4(x, t = 0) = 0. This envelope becomes popu- 
lated as a result of the coherent 4WM process. Substitut- 
ing the SVEA form for the wavefunction into the GPE, 
collecting terms with the same phase factors, multiplying 
by the complex conjugate of the appropriate phase fac- 
tors, and neglecting terms that are not phase matched 
(those for which momentum and energy are not con- 
served) we obtain a set of coupled SVEA equations for 
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Only phase-matched terms, for which kj 
k r — k s = and ujj qrs = loj — uj q + u> r — ui s = 0, are 
retained on the right hand side of Eqs. (||). 

For simplicity, we consider explicitly the SVEA equa- 
tions for the case where only two central momentum com- 
ponents, and 2Kkph, are present. Then, only "phase- 
modulation" nonlinear self- and cross-energy interaction 
terms are present, as opposed to the case when three cen- 
tral momentum components are present and 4WM terms 
also arise. It is convenient to go to a center of mass frame 
moving with velocity v cm = Kk p h/m (see Fig. |l|). In this 
frame, wavepacket 1 has momentum —Kk p h, wavepacket 
2 has momentum Kk p h, and elastic scattering from the 
±?ikph wavepackets creates a spherical shell expanding 
in 47r steradians with momentum |?ik e ; as | = |?ikp^|. The 
SVEA equations in this frame are given explicitly by: 
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Eq. (|4|) explicitly separates out the fast oscillating phase 
factors representing central momentum Kkj and kinetic 
energy Ej — fitoj — fik 2 /2m and defines the slowly 
varying envelopes <&j, which vary in time and space on 



where a = 8iraQ and v re i — 2Tik p } l /m = 2v cm is 
the relative velocity of the two wavepackets. The fac- 
tor of 2 multiplying the nonlinear cross-energy interac- 
tion terms results from expanding |\l/| 2 \Ij with ^(x, t) = 
^2j = i ^(x, t) exp(ikj-x — iujjt) and collecting the phase- 
matched terms appropriately. In the self-energy term 
proportional to |$i| 2 $i, i = 1,2, only the lead term 
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in Eq. (@) is retained. However, both terms in Eq. (g) 
are retained in the cross-energy terms |$j| 2 $j, i = 1,2, 
j = 2, 1, leading to the elastic scattering loss terms pro- 
portional tO V re l<7. 

The form of the elastic collisional loss terms can also 
be motivated by a classical hydrodynamic picture of elas- 
tic collisions of a cloud of atoms having a central velocity 
vi = Pi/m and density ni(x, t), with a cloud of atoms 
having velocity v 2 and density ri2(x,t). The atomic den- 
sities can be determined from: 



<9ni(x, t) 
dt 



vi • Vni(x,t) = -|vt - v 2 |<7mn 2 , (8) 



and a similar equation for a " 2 gj* ■ Thus, elastic scat- 
tering takes atoms out of both clouds. For sufficiently 
slow relative atomic velocities so that only s-wave scat- 
tering occurs, and both atoms have the same spin quan- 
tum numbers, a = 87retQ. Here the relative velocity is 
v re i = |vi — V2I = 2v cm . The extra factor of | in the 
loss terms in Eqs. are due to the fact that these 

are equations for amplitudes ($), not densities (l*^ 2 )- 
Inelastic scattering with cross section crj„, if present, is 
easy to include by replacing a by a + a% n . 

If we apply this theory to condensates with two spin 
components, the cross section which appears in the loss 
term due to collisions between wavepackets of the two 
different components is a — 47raQ, as expected from two- 
body scattering theory for different spin components. 
In the SVEA derivation this follows from the fact that 
the cross phase modulation terms are then of the form 
-±U N T (\m 2 )^, rather than -j-U Q N T (2\^\ 2 )^ as 
here (note the factor of 2), where the superscripts a and 
b denote spin indices. 

As a first example, we consider a condensate of Na 
atoms in the F = 1, M = — 1 Zccman sublevel in 
a cigar-shaped trap elongated in the z-direction. Af- 
ter the trapping potential is turned off, the condensate 
is allowed to freely evolve for 600 /is, and a short du- 
ration Bragg scattering pulse is applied that creates a 
27ikp/j = (2h/\ p h)z momentum component. We con- 
sider the case with half the initial atoms in the high 
momentum component and half in the parent conden- 
sate. We evolve the condensate wavepackets using a full 
3D implementation of Eqs. (||)-(0) until the wavepack- 
ets move apart and are physically separated. Fig. (||) 
shows Nj/Nt versus the aspect ratio, R aS pect, of the ini- 
tial elliptically shaped BEC for two different initial total 
number of atoms in the BEC, Nt — 1.0 x 10 6 atoms 
and Nt — 3.0 x 10 6 atoms respectively. Here Nj is 
the total number of atoms remaining in both the Ohkph 
and 2h\iph wavepackets after the wavepackets separate. 
Thus, Nj/Nt = 1 — L where L is the fractional loss of 
atoms from the mean-field due to elastic scattering. The 
Thomas- Fermi aspect ratio is related to the trap frequen- 
cies by Raspect = xtf/ztf = w z /lj x . The actual trap 
frequencies in our calculation were v z = 30.7 Hz, and 
v x = v y = v z /Raspect- The figure shows that the loss 



increases as the aspect ratio decreases, and as the total 
number of atoms increases, reaching 60% for 3 million 
atoms and R asp ect ~ 1/20. Nj/Nt rises slowly to unity 
as Rasped gets large. 
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FIG. 2. Fraction of atoms remaining in the condensate 
wavepackets after the wavepackets have separated. A fast 
daughter wavepacket with half the initial number of atoms 
and with momentum 2hk in the z-direction moves through 
the remaining parent condensate with zero central momen- 
tum. The dashed line shows the result of the heuristic model 
for 10 6 atoms. 



A simple heuristic model helps to explain the mag- 
nitude of the losses. Assume a uniform atom density 
of 71 = Nt/V in a cylinder of length I — 2ztf, area 
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and volume V = IA. If an equal num- 
ber of atoms Nt/2 is assumed to be in the daughter 
and parent wavepackets, and depletion of n during the 
interaction is ignored, a simple argument shows that 
the fraction of atoms remaining after the packets sep- 
arate is (£ m j p /£)(l — e~ l / lm fv), where £ m j p = (f ct) -1 is 
the mean-free-path for the collision. For example, with 
v z = 30.7 Hz, Rasped =0.1, and Nt = 10 6 atoms, we 
find that I = 120 /im and l m j v = 140 /xm are comparable 
in magnitude. The dashed line in Fig. |^ shows that this 
simple model qualitatively accounts for our results. 

In the NIST 4WM experiment g, the Na(F = 1, M = 
— 1) condensate is exposed to Raman scattering pulses 
which create copies of the parent condensate at central 
momenta Kk.2 = (h/\ p h)(% + y) and Kk^ — (2h/\ p h)%, 
leaving part of the atoms in the Tiki = wavepacket. 
The treatment of elastic scattering from the disparate 
momentum components of the wavepacket in the 4WM 
experiment is similar to the description above for the two 
momentum component case. Now, there are three elastic 
scattering loss terms for each SVE momentum compo- 
nent $j arising due to the cross-phase modulation terms 
of each momentum component with the other three mo- 
mentum components. We also included the momentum 
dependent correction term in Eq. (^|) in the coupling con- 
stant for the 4WM source terms on the right hand side of 
Eq. (||); it only slightly decreases 4WM at large Nt- In 
the experiment a trap with v x = 84 Hz, v y = 59 Hz and 
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Vg = 42 Hz contained a Na BEC without a discernible 
non-condensed fraction. The condensate was exposed to 
Raman scattering pulses 600 /j,s after the magnetic har- 
monic potential was turned off. Fig. || shows the fraction 
of atoms in the 4WM output wavepacket as a function 
of the initial total number of atoms Nt as determined 
(1) experimentally (circles), (2) by calculation without 
including elastic scattering loss for a ratio of atoms in 
the three initial wavepackets of 7:3:7, and (3) by calcu- 
lation including elastic scattering. The effects of elastic 
scattering are pronounced for large values of Nt, with 
the percent loss due to elastic scattering reaching 44% 
for 5 million atoms. The discrepancy with experiment is 
reduced significantly by including loss due to elastic scat- 
tering. Possible remaining sources of discrepancy include 
micromotion of the BEC in the time-orbiting-trap, laser 
misalignment and a small finite temperature component 
of the BEC. 
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FIG. 3. Fraction of atoms in the 4WM output wavepacket, 
Na/Nt, versus the number of initial atoms, Nt- The dots 
represent experimental data [3, the solid curve is the 3D cal- 
culation without elastic scattering, and the heavy solid curve 
is the 3D calculation with elastic scattering. 



It is useful to put the elastic scattering discussed here 
into perspective. The mean-field wavefunction for the 
zero temperature BEC is a symmetric product of identi- 
cal orbitals, each orbital being a coherent superposition 
of momentum wavepackets. Elastic scattering between 
the various momentum components of a BEC results in 
atoms which can not be described by the mean-field since 
the "modes" into which the atoms are scattered (there 
are an infinite number of scattering angles, or modes, to 
scatter into) are not macroscopically populated p| , |l7j |. 
The momentum components in the spherical shell can 
not be generated using the dynamics of the GPE or the 
SVEA equations since neither contains terms that pro- 
duce such momentum components. However, the SVEA 
equations do allow collision losses to be treated. In con- 
trast to the formation of the spherical shell of elastically 
scattered atoms, the fourth wave in 4WM is explicitly 
generated by the GPE or the SVEA equations. 

The scattering of atoms into the spherical shell and 
the loss of atoms from the condensate are a result of 



Hamiltonian dynamics; no interactions with a bath, and 
therefore no incoherent processes described by T± or T2 
relaxation times, are necessary. Our treatment of the 
process has, for convenience, used an imaginary poten- 
tial that serves as a mechanism to take atoms that are 
elastically scattered out of the condensate. It must also 
be noted we were able to carry out our procedure for 
modeling the loss of atoms from the condensate only as 
a result of making the slowly varying envelope approxi- 
mation, by which we could track the density of atoms in 
each momentum wavepacket individually. As the elastic 
scattering loss increases, further scattering of the elas- 
tically scattered atoms with the condensate atoms will 
become increasingly important. This mechanism is not 
included in our treatment and would require following 
the dynamics of the elastically scattered atoms in detail. 
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